Mode Coupling in Quantized High Quality Films 



Yiying Cheng and A. E. Meyerovich 
Department of Physics, University of Rhode Island, 
2 Lippitt Rd., Kingston 
RI 02881-0817 

The effect of coupling of quantized modes on transport and localization in ultrathin films with 
quantum size effect (QSE) is discussed. The emphasis is on comparison of films with Gaussian, 
exponential, and power-law long-range behavior of the correlation function of surface, thickness, or 
bulk fluctuations. For small-size inhomogeneities, the mode coupling is the same for inhomogeneities 
of all types and the transport coefficients behave in the same way. The mode coupling becomes 
extremely sensitive to the correlators for large-size inhomogeneities leading to the drastically distinct 
behavior of the transport coefficients. In high-quality films there is a noticeable difference between 
the QSE patterns for films with bulk and surface inhomogeneities which explains why the recently 
predicted new type of QSE with large oscillations of the transport coefficients can be observed 
mostly in films with surface-driven relaxation. In such films with surface-dominated scattering the 
higher modes contribute to the transport only as a result of opening of the corresponding mode 
coupling channels and appear one by one. Mode coupling also explains a much higher transport 
contribution from the higher modes than it is commonly believed. Possible correlations between 
the inhomogeneities from the opposite walls provide, because of their oscillating response to the 
mode quantum numbers, a unique insight into the mode coupling. The presence of inhomogeneities 
of several sizes leads not to a mechanical mixture of QSE patterns, but to the overall shifting and 
smoothing of the oscillations. The results can lead to new, non-destructive ways of analysis of the 
buried interfaces and to study of inhomogeneities on the scales which are inaccessible for scanning 
techniques. 

I. INTRODUCTION 

Progress in material technology, especially in nanofabrication, ultrathin film deposition, ultraclcan and high vacuum 
systems, etc., requires better understanding of the effect of remaining bulk inhomogeneities or surface defects on 
physical processes in high-quality systems. In high-quality systems, these remaining inhomogeneities small and smooth 
with a low amplitude and a relatively large lateral scale. In some cases, such as in ultrathin films, the lateral scale 
of bulk and surface inhomogeneities can even be much larger than the film thickness. Scattering by such small, but 
long-range inhomogeneities is crucial for transport in ultrathin and/or clean systems in which the particle mean free 
path is comparable to the system size. 

Below we consider the effect of random, mostly large-scale, bulk, surface, and thickness fluctuations on quantum 
transport in quantized quasi-2Z? systems such as quantized flow channels, waveguides, or ultrathin metal films. We 
will look at the single-particle diffusion coefficient D in a channel as a function of the particle energy and the channel 
width and at the low-temperature mobility fi (conductivity a) as a function of the film thickness and the Fermi 
wavelength. The main issue is to find how sensitive is the transport to the statistical properties of inhomogeneities, 
i.e., to the structural or thickness fluctuations with small amplitude and large correlation radius. Here we have in 
mind large-size surface steps and thickness fluctuations for ultrathin films, slow long-range bending of fibers or films, 
slowly fluctuating bulk fields, etc. One of the main goals is to separate the effect of the scattering-driven mode 
coupling from other scattering effects. 

The choice of quasi-2£> systems is explained by a desire to avoid divergence of surface fluctuations and strong 
localization effects which are inherent to ID systems. In contrast to II? systems, the randomly fluctuating 2D 
: surfaces are stable while the localization length in systems with weak fluctuations is exponentially large. 

Usual approaches to bulk and surface fluctuations are different from each other. The bulk fluctuations are routinely 
described via the fluctuating bulk potential V (r) or, whenever possible, via the scattering T-matrix, T (p, p'). Since 
V and T are tied to each other via the integral equation, 

f = V + TGV (1) 

(G is the Green's function), these two descriptions arc, in principle, equivalent (and, in the case of weak fluctuations 
- identical). Below we assume that the bulk inhomogeneities are defined by their scattering T-matrix and that this 
T-matrix is known. 

The prevalent way to characterize the surface roughness or thickness fluctuations is to use the correlation function 
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of the surface inhomogeneities 

C (s) = C (|s|) = (e(si)C(si + s)> = A- 1 J e(si)C(si + s)ds U (2) 

where s gives the 2D coordinates along the surface, £(s) describes the deviation of the position of the surface in the 
point with 2D coordinates s from its average position, (£(s)) = 0, and A is the averaging area. This equation assumes 
that the correlation properties of the surface do not depend on the lateral direction. 

Any transport theory should yield an explicit dependence of the transport coefficients (the particle mean free path, 
diffusion or mobility coefficients, lateral conductivity, etc.) on the correlator of surface inhomogeneities Q (s) and/or 
the bulk scattering T-matrix. Both of these objects can be characterized by their average amplitudes £, types of 
decay (exponential, power law, etc.) at large distances or momenta, and the correlation radii R which characterizes 
this decay. For the surface correlation functions £, parameters I and R are the average "height" and "lateral size" 
of surface inhomogeneities or thickness fluctuations. For the bulk scattering T-matrix, I is the average amplitude 
in the dependence of the scattering amplitude on dimensionless momentum pR/h and R provides a scale for decay 
of the scattering amplitude at large momentum transfers. The transport coefficients are determined by the relation 
between the particle wavelength A, the width of the channel L, and the correlation radius of inhomogeneities R. If the 
fluctuations are weak, the fourth length parameter, £, (more precisely, its square) enters as a coefficient. For example, 
the conductivity a of degenerate fermions and the single-particle diffusion coefficient D can be parameterized as 

2e 2 R 2 

° = - T ^-/(A,L,J2) ) (3) 

D = ^ d (A,L,R). (4) 

with virtually identical functions / and d. The reason for this commonality is that both of these transport coefficients 
are expressed via the same combination of the zeroth and first angular harmonics of the scattering probabilities. 

We consider ultrathin systems with quantized motion across the film. Scattering by surface and bulk inhomo- 
geneities could cause coupling of the otherwise distinct quantum modes. We will show that the qualitative behavior 
of the transport coefficients is extremely sensitive to this scattering-driven coupling of the modes which, in turn, is de- 
termined by the long-range behavior of the structural or thickness fluctuations. Mode coupling has already attracted 
considerable attention for localization lj and transport, especially in the context of applications of the random matrix 
theory Though most of the applications involved transport in systems with bulk disorder, systems with surface 
disorder were also considered 

What has been mostly ignored is the sensitivity of the mode coupling and, in the end, transport, to the type 
of correlation behavior of disorder. It has been often assumed that the correlation function is short-range (5-type 
correlations, hard spheres, etc.) without long tails. What is more, in the case of impurities a natural assumption is 
that the correlation radius R ("size"} of the disorder is relatively small. Under these assumptions the mode coupling 
is indeed featureless (though robust) and does not lead to any striking effects which depend on the nature of disorder. 
However, in high-quality quasi-2D samples, it is possible to observe bulk and surface disorder with various correlation 
properties and fluctuations of different scales. In this case, as we will see below, the mode coupling and transport can 
follow several distinct scenarios. These different types of behavior are determined by the rate of decay of correlations 
and their scale and not by the origin of fluctuations or nature and spectra of the particles and waves. 

Recently we predicted |5( a new type of quantum size effect (QSE) with huge, large-period oscillations of conductivity 
a (L) in films with dominant surface scattering. This surface-driven effect is so large that a real puzzle is why this 
new type of QSE has not been observed earlier in high quality films with bulk scattering. Below we will answer this 
question by comparing films with bulk and surface scattering. We will also analyze the contribution of different modes 
to transport. A usual assumption is that in films with surface scattering the main contribution to transport comes 
form the grazing particles - modes with the lowest quantum numbers. We will see that the scattering-driven mode 
coupling makes this assumption wrong even for the highest quality films. The last important issue that will be studied 
below is the possibility of extracting information on the quality of the film surfaces, including the buried ones, from 
the non-destructive transport measurements. Interestingly, QSE in transport can even reveal a possible correlation 
between random surface inhomogeneities from the opposite walls. 



II. TRANSPORT IN QUANTIZED SYSTEMS 



Below we deal mostly with the conductivity of metal films and the single-particle diffusion in quasi- 2D channels. 
QSE experiments in metal films include conductivity [Tol Il6|. spectroscopy HJ, susceptibility O, and STM [l3| 
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measurements. The signature feature of QSE in metals is a pronounced saw-tooth dependence of the lateral con- 
ductivity on, for example, film thickness, cr(L), common for both bulk [l4j and surface [l5j scattering. However, 
the QSE experiments in metals have to overcome a difficulty which one does not encounter in semiconductors. The 
period of the saw-tooth QSE oscillations in the dependence a (L) is usually small, nearly atomic, h/pp, making the 
saw-tooth behavior of transport almost impossible to observe. For this reason, typical experimental objects are lead 
or semimetal films such as bismuth. In contrast to this "standard" saw-tooth effect, QSE, which is inherent to high 
quality films [f|, results in smooth, large-period oscillations of a (L) at relatively large values of pfL/H. This QSE 
could be observed in a wider group of metals. Large-period QSE oscillations have already been observed (see the 
second Ref. [13); however, the experimental details are still sketchy. This issue is also related to the long-standing 
controversy on the influence of the structure of the nanoscale film on its resistivity |l6j| . 

Another group of seemingly different physical applications involves the single-particle diffusion in quantized flow 
channels. Though the typical examples - quantized helium quasiparticles in ultrathin channels [TtI Il8( and ultracold 
neutrons in gravitational traps |l9j| - seem to be far apart from metal films, the descriptions of the transport processes 
in such diverse systems are very similar to each other. 

Recently, we developed a transparent semi-analytical formalism for transport in systems with rough boundaries that 
allows simple uniform calculations in a wide range of systems and for various types of roughness |J, |20j . One attractive 
feature is that this formalism treats the scattering by surface inhomogeneities using the same transport e quat ion as 
for the bulk imperfections and even allows one to study the interference between bulk and surface scattering |2l| . This 
formalism unites earlier approaches by Tesanovic et al [23, Fishman and Calecki [23|. Kawabata |24|. Meyerovich 
and S. Stepaniants [25|. and Makarov et al [^3 (for a brief review of different theoretical approaches see Refs.Hlll). 
In this paper we apply this approach to the study of the fluctuation-driven coupling of quantized modes. The limits 
of applicability of our approach are discussed in detail in Refs. 0, [^J . 

Since the 2D mobility of particles is described by essentially the same equations as the exponent in the expression 
for the localization length in films, our study also provides the dependence of the localization length on the type of 
the correlation function of random surface and bulk inhomogeneities. 

The paper has the following structure. In the next Section we introduce the transport equation and expressions 
for the transport coefficients. Section IV briefly describes various types of the surface inhomogeneities and bulk 
T-matrices used in the computations. The results are given in Section V and conclusions - in Section VI. 



III. TRANSPORT EQUATION AND TRANSPORT COEFFICIENTS 



QSE is caused by quantization of motion in the direction perpendicular to the film, p x — > njh/ L, which leads to a 
split of the energy spectrum e (p) into a set of minibands, e (p x , q) — > e (irjh/L, q) = Cj (q). For simplicity, we consider 
mostly particles with a parabolic spectrum, 



£ j (q) 



i 

2m 



(irjh/L) 2 + q 2 



(5) 



The extension to non-parabolic spectra is discussed in Section VE. We will look at two similar transport problems, 
namely, at conductivity of degenerate fermions, 



qj = qFj 



2mep — (irjh/Ly 



1/2 



(6) 



{ep = ir 2 h 2 /2mKp is the overall Fermi energy, is the Fermi wavelength, and qj (ep,L) is the Fermi momentum 
for the miniband j) and at single-particle diffusion for particles with energy E = ir 2 h 2 /2mA 2 , 



1j 



2mE - {■Kjh/Lf 



1/2 



(7) 



where qj (E,L) is the momentum of such a particle in the miniband j. Both problems are computationally similar 
and reduce to almost identical sets of linear equations . 

In the case of conductivity of degenerate fermions, the transport equation for the distribution functions rij (q) , 



d 2 q' 



(2irhy 



(8) 
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reduces, after standard transformations, to a set of linear equations 



qj/m = v r {q r ) /T jf , (9) 

f 

2 



T "' i" 

where = Vj5{e — ep)eE is the first angular harmonic of the distribution function rij (q) at q = qj, and 

W^/' 1 ' (QjjQf) are the zeroth and first harmonics of W (qj— q^) over the angle q^q^/- The solution of Eqs. 
provides the conductivity of the film: 

e 2 x - 

a = ~^22 v ^ (<lj)qj- (10) 

3 

The matrix r is diagonal when the scattering-driven mode coupling is negligible with respect to the in-band scattering. 
The single-particle diffusion is similar to the conductivity problem for degenerate fermions. In equilibrium, a particle 

with energy E can be in any of S accessible minibands 6j (q) = (l/2m) (irjh/L) 2 + q 2 for which ej (q = 0) < E. 
The equilibrium distribution function n(°> is 



,.(o) 



and the transport equation reduces to a set of S coupled linear equations for distributions Uj (qj) with momenta q 3 
which is almost identical to Eqs. © , ifTU|) : 

j^lj^P = -J^i'tejOAtt'. ( 12 ) 
i' 

1 5 

3=1 

where Vp is the density gradient that causes the diffusion, I? is the single-particle diffusion coefficient, and r is still 
defined by Eq.Q. The single-particle mobility coefficient b is related to D by the Einstein equation D = bE. 

The results for the single-particle diffusion also provide the mean free path C and the exponent in the expression 
for the localization length 1Z that describes localization caused by particle scattering by random wall and/or bulk 
inhomogeneities 

1l = Cexp[nmSD/h]. (13) 



IV. SURFACE CORRELATION FUNCTION VS. BULK SCATTERING AMPLITUDE 

In the case of surface or thickness fluctuations, we consider an infinite 2D channel (or film) of the average thickness 
L with random rough boundaries 

x = L/2- ^(y, z), x = -L/2 + £ a (y, z). (14) 

(the walls are assumed hard with infinite potential). The inhomogeneities are small, £j 2 {Vi z ) ^ L, and random with 
zero average, (^) = (£ 2 ) = 0. Their correlation function Q ik (s) and its Fourier image £ ifc (q), which is often called 
the power spectral density function or, in short, the power spectrum, are defined as 

(N) = fe(si)^(si + s)) = A- 1 J &(si)£ fc (si + s)d Sl , (15) 
C«(|q|) = / ' d 2 S e^X tk (\s\) = 2* Qk (s) J Q (qs) sds 

where s = (y, z) and q = (q y , q z ) are the 2D vectors. In homogeneous systems, the correlation function depends only 
on the distance between points |si — s 2 and not on coordinates themselves. The correlation functions £ n and £ 22 
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describe intrawall correlations of inhomogeneities, and £12 — C21 " the interwall correlations. Usually, but not always, 
the inhomogeneities on different walls are not correlated with each other, £ 12 = 0. Thus, everywhere, except for 
Section VF, it is assumed that £12 = 0. To avoid parameter clutter, we also assume that the correlation parameters 
are the same on both walls, Cn = C22 = C- Then the effective correlator contains 2£ (s) with £ (s) given by equations 
below. 

Surface inhomogeneities exhibit a wide variety of types of the correlation functions ■ To have a meaning- 
ful comparison, we consider the correlation functions that involve only two characteristic parameters, namely, the 
amplitude (average height) I and the correlation radius (average size) R of surface inhomogeneities. 

The most commonly used in theoretical applications is the Gaussian correlation function, 

C (s) = i 2 exp (s 2 /2R 2 ) , C (5) = 2ir£ 2 R 2 exp (-q 2 R 2 /2) (16) 

including its limit for small correlation radius R — > 0, i.e., the <5-type correlations, 

C (s) = £ 2 R 2 S (s) /a, C (g) = 2tt£ 2 R 2 . (17) 

Sometimes, a better fit to experimental data on surface scattering is provided by the use of the exponential corre- 
lation function 

2tt/ 2 R 2 

C (-) = e exp (s/R) , C («) = —3^ (18) 

(1 + q z R z ) ' 

or by the even more long-range, power-law correlators 

C ( s ) = 2 ^ C ( q ) = 2tt£ 2 R 2 {qR) ^ (qR) (19) 

sw (l + s 2 /i? 2 ) 1+A1 2»- 1 T(fj) ' y ' 

with different values of the parameter /i. The asymptotic behavior of the power spectra, i.e., of the functions (qR), 
is essentially exponential. The most commonly used are the Staras function with /i = 1 and the correlator with fi = 1/2 
which has the plain exponential power spectrum £ (g) , 

C (q) = 2t:£ 2 R 2 exp (-qR) . (20) 

The last class of correlation functions covers the power-law correlators in momentum space, 

C f q ) = 27r£2R \^ , ( (a) = £ 2 } S/R)X , K X (a/R) , (21) 

(l + g 2 i? 2 ) 1+A 2 A r(l + A) V ' V ' 

which are exponential functions in the configuration space. The correlators from this group include the Lorentzian in 
momentum space A = that was observed in Ref.0 and the exponential correlator l|18|) at A = 1/2. 

The constants in all these correlators are chosen in such a way that the value of ( (q = 0) = 2n£ 2 R 2 is the same. 
This provides a reasonable basis of comparison for transport coefficients in films with all these different types of 
random surfaces. Indeed, the scattering cross-section for q — > does not depend on the details of short- and mid- 
range structure of surface inhomogeneities. Therefore, at Fermi momenta qF — > (more precisely, at qpR <SC 1), the 
transport coefficients should be the same for all random surfaces. (The only exception is the Lorentzian (??) for which 
£ (q) diverges at small q; however, the Lorentzian is non-physical anyway). Some useful analytical expressions for the 
angular harmonics of these correlation functions can be found in Ref . [5j . 

To have a uniform comparison of the results throughout the paper, we will plot numerical results for the Gaussian 
correlator, the power-law correlator in configuration space l|19[) with fi = 1 (the Staras function), and the power-law 
correlator in momentum space (|211 with A = 0.5 (exponential correlator in real space). 

In the case of bulk modulation, it makes sense to start directly from the bulk scattering amplitude T (p, p') which, 
in the case of quantized films, is transformed into the matrix 

T (p, p') = T (p x , q;p' x , q) -» T jf (q, q') = T (q 3 - q' f ) . (22) 

For this matrix, we also looked at the Gaussian form similar to (I16fl . power law form similar to (|21|l with A = 0.5 
(standard Debye screening with an exponent in real space) , and exponential form l|19l) with ^ = 1 . Another interesting 
possibility here would be an oscillating function which corresponds to an alternative type of screening in metals. 
However, such a function would require us to introduce two lateral length parameters which would make a meaningful 
comparison between the correlation functions impossible. 
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In what follows we compare the transport properties of the films <|16|) - HP in a wide range of film thickness L, 
correlation radius R, and particle wavelength A = n/q (or the 2D particle density N). 

The transport equation for both bulk and surface imperfections is the same, Eqs. (JHJ - (|12|l The only distinguishing 
feature is the dependence of the scattering probabilities Wjf (q, q') on the correlation function of surface fluctuations 
C(|q-q'l) EH. 



Wjf (q,q') = 



>-L 2 I 



Cn +C22 + 2Ci 2 



1L 

L 



and the T-matrix for the bulk imperfections, 

Wjf (q,q') 



(23) 



(24) 



Note that in contrast to Eq. (|23|l . the dependence of W (I24|l on band indices j,f is generally not known explicitly 
and is determined by the T-matrix (|22f) . This dependence is the same as in Eq. l|2.'ill when the fluctuating bulk field 
is essentially two-dimensional and can be factorized as 



U (x) + (dU/dx) £ (s) 



(25) 



where ians are the coordinate across and along the film and U (x) is the field without fluctuations. This situation, 
though realistic, is by no means general. If, for example, the bulk field fluctuates only along the film and does not 
change across the film, then T [qj — q^,) = SjfT (q — q'). The off-diagonal elements of T (q 3 — q^/) and, therefore, 
the mode coupling are associated only with the variation of the bulk field across the film. All this makes modeling 
of the function T (q^ — q^) more ambiguous than for the surface correlator. There are two ways of dealing with this. 
The first one is to assume that the bulk fluctuations in ultrathin films have the form l|25|) and, essentially, to use 
the expression similar to Eq. (|23(l for the scattering probabilities (of course, without the oscillating interwall term 

C12 )■ Another approach, which is more appropriate for thicker films is to assume that the bulk fluctuations 

are truly three-dimensional and are not affected by the film boundaries except from the straightforward quantization, 



T r y (q, q') = T (irjH/L, q; nj'h/L, q') . 
Then, for example, the Gaussian fluctuations in momentum space lead to the following form of W: 



„r , n 8ir 5 h 5 £ 2 R 2 , 2 ~ ,^ 
Wjf (q, q') = m2£6 exp {-q 2 R 2 /2) exp 



2L 2 



(26) 



(27) 



where we introduced parameters £ and R in such a way so that to make the description as close as possible to the one 
with the Gaussian thickness fluctuations (ftT))l (or with the one with the Gaussian bulk fluctuations of the type 



Wjf (q,q') 



8n 5 h 5 l 2 R 2 
m 2 L e 



j 2 j' 2 exp(-q 2 R 2 /2 



(28) 



V. RESULTS AND DISCUSSION 



General comments 



As it is mentioned in Introduction, the 2D conductivity a of the film has the dimensionality of conductance, 
e 2 /?i.The conductivity depends on the relation between three length scales - particle (Fermi) wavelength A^ = 1/qF, 
the width of the channel L, and the "lateral size" (correlation radius) of inhomogeneities R. The fourth length 
parameter, £, is perturbative and enters conductivity as a perturbative coefficient l/£ 2 , 



2e 2 R 2 



f(L/A F) R/L) 



(29) 



Note, that the conductivity diverges in the limit of vanishing inhomogeneities £ — > or R 
The single-particle diffusion coefficient D can be parameterized in a similar way, 



D 



f^R 2 

m £ 2 



d(L/A,R/L) 



(30) 
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where the dimensionless function d is given by the same combination of the zeroth and first harmonics of the scattering 
probability W, 

d(L/A,R/L) = ^f(L/A,R/L) (31) 

and S = Int(L/A) is the number of quantized energy minibands accessible to the particle with the wavelength A. The 
presence of the step- wise function S(L/A) in the denominator of Ea. l|31(l can lead to a visible difference in shapes 
between functions d (L/A, R/L) and / (L/A, R/L): at large R/L the function / is practically smooth while the function 
d has residual singularities in the point where the number of accessible minibands changes. These singularities can 
serve as useful markers that help to separate the "standard" QSE and the QSE of Ref. ||. Therefore, we will plot the 
results for d (L/A, R/L) more often than for / (L/A). Apart from these small-amplitude singularities, the functions 
d and / are similar. 

Below we will plot the dimensionless functions d (L/A) and / (L/Ap) at various constant values of R/L. The orders 
of magnitude of these functions change widely depending on the type of surface and bulk correlators and the value 
of R/L. In order to provide a better visual comparison between these functions, we will normalize d (x) and / (x) by 
their (usually, maximal) values at the maximal value of x — L/A in the calculation. In other words, we will plot the 
normalized functions 

d(x) /d(x max ) , f(x)/f(x max ) (32) 

where the values of the coefficients d (x m&x ) and / (a; max ) are irrelevant for our purposes. In this way, all the functions 
change from at x = to 1 at x = £ max and can be visually compared with each other. The values of x max in our 
computations vary from 35 to 100 which means that there arc between 35 and 100 quantized minibands accessible to 
the particles. 

The data below involve the Gaussian correlator (|16fl . power-law correlator (|19fl with // = 1 (the Staras function), 
and the exponential correlator in real space l|18f) (i.e., the correlator l|21|) with A = 0.5). We will refer to the latter 
two as the power-law and exponential correlators without specifying /i and A. 



B. Mode coupling 

To underscore the effects of mode coupling, we start from the calculation with the artificially frozen mode coupling 
(the off-diagonal components of the matrix of scattering probabilities Wjj> are not calculated, but made equal to 
zero, Wjji = WjSjj'). This is a good starting point since in high quality films with R/L 3> 1 the mode coupling is 
suppressed anyway. In this case the transport equations © and i|12|) can be solved analytically [2(i|: 



2e 2 ^ q 2 



3h 2 m 2 - U/W 

and similarly for diffusion, 



(33) 



2 

J ^4^E (J' (i) ' ( 34 ) 

where W^ ' 1 ^ are the zeroth and first angular harmonics of the transition probabilities 1231) , l)24[l. Note, that since 
Wj for surface correlators grows proportionally to j 4 , Eq. (|23|l . the sum in Eq. (|33J) is rapidly convergent and 
the contribution of the higher modes in systems without mode coupling can be negligible. This means that for the 
surface scattering without mode coupling the dependence of the conductivity a on the film thickness is practically 
a smooth function, Figure 1 (small kinks on the curves cannot be seen with the resolution of the figure). This is 
not so for the single-particle diffusion (|34l) which contains a step- wise factor S(L), Figure 2. The small saw-tooth 
drops on the curves d(L) at the points in which S = Int (L/A) changes by 1 (quantum size effect, QSE) can serve as 
useful markers on the curves that help to separate the standard QSE from other phenomena. Therefore, more often 
than not we will plot d (L) rather than a (L). Figures 1 and 2 combine the data for the Gaussian, exponential, and 
power-law inhomogeneities and five different values of R/L = 0.1; 1; 10; 50; 100. All normalized curves / (x) / f (50) 
and d(x) /d(50) are identical to each other except, of course, for the normalization coefficients / (50) and d (50), 
which change by the orders of magnitude depending on R/L and the type of inhomogeneities. 
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C. Bulk vs. surface scattering 

In Ref. we reported the existence of a new class of QSE for a "boring" type of high-quality films R ^s> L 
with Gaussian or exponential (in momentum space) surface and thickness fluctuations. This QSE manifests itself as 
giant oscillations of a (L) with a relatively large period which is directly related to the correlation radius of surface 
inhomogeneities, 

L S ~ |V(2j + l)i?A, (35) 

where Lj are the positions of the peaks in a (L). 

The effect is so large and well pronounced that a natural question is why has not it been observed earlier in, for 
example, impurity scattering. As we discussed in the Introduction and Section IV, the equations that describe surface 
and bulk scattering are so similar that it looks like this new QSE should be observed in bulk scattering as well. 
Surprisingly, this is not the case. Figure 3 presents the normalized diffusion coefficient for surface and bulk scattering. 
Both surface and bulk correlation functions are Gaussian with the same large value of R/L = 50. The curve with 
surface scattering exhibits large QSE oscillations with peaks in positions l|35|). while the bulk scattering results in 
a monotonic dependence similar to that for systems with frozen mode coupling in Figures 1 and 2. This surprising 
difference between the bulk and boundary scattering requires an explanation. 

Our explanation of the "new" QSE, Eq. , in Ref. @ was that in high-quality films with R/L ^> 1 the mode 
coupling is largely suppressed because the scattering-driven change in momentum Sq ~ l/R is insufficient to induce 
transitions between the modes which are widely separated between themselves at small L, Se oc 1/L 2 . The mode 
coupling processes turn on one by one only at the values of thickness given by Eq. i|35|) . This should be the same for 
both surface and bulk scattering. What is not the same is the effect that this opening of the mode coupling channels 
has on transport. In the case of surface scattering, the modes with the lowest quantum numbers (grazing particles) 
make the largest contribution to transport (without the quantum cut-off, the contribution from the particles in the 
lowest mode - grazing particles moving parallel to the surface - would have been infinite). Therefore, the threshold 
opening of individual mode coupling channels for the lowest modes, which almost doubles the corresponding cross- 
sections, is very noticeable in transport. In the case of bulk scattering, the situation is different. Here all modes 
contribute more or less equally to transport. Therefore, in the case of a large number of available modes, the opening 
of few mode-coupling channels in points l|35l) is unnoticeable and the transport coefficients remain nearly the same as 
in the pure diagonal case. This explains why this new QSE (I35|l has not been discovered earlier when studying the 
bulk-dominated transport. 

Figure 4 illustrates very different sensitivity of the surface- and bulk-driven transport to mode coupling. The figure 
contains the normalized single-particle diffusion coefficient in the cases of bulk and surface scattering. In both cases 
the inhomogeneities are Gaussian with four different values R/L — 0.1; 1; 10; 50. In the case of surface scattering, the 
pattern evolves from the "usual" QSE saw-tooth curve for R/L = 0.1 to the new QSE with huge oscillations on more 
or less smooth curve for R/L = 50 (the curves are marked by the value of R/L). In the case of bulk scattering, all 
the curves split into two groups of smooth curves for small and large R/L (with, correspondingly, robust and mostly 
suppressed mode coupling). Though the shapes of the curves from these two groups are distinctly different, the 
difference in shapes, in contrast to surface scattering, is rather quantitative than qualitative. Of course, the saw-tooth 
anomalies, which are inherent to QSE, are more distinct on the curves with robust mode coupling at small R/L. 

Summarizing, QSE and the manifestations of the mode coupling are distinctly different in surface- and bulk dom- 
inated quantized systems. This is due not to the difference in mode coupling processes themselves, but due to the 
different patterns for contributions from individual modes, especially for the grazing particles. This also answers a 
puzzling question why the new type of QSE is observed primarily in systems with the dominant surface relaxation. 

D. Opening of mode coupling channels and mode contributions for surface scattering 

The next step is the analysis of contributions from individual modes and from the mode coupling transitions. 
We start from the data for the artificially frozen mode coupling transitions in Figures 1,2 for surface scattering 
and turn on such transitions. The results strongly depend on the size of inhomogeneities R/L. Analysis of the 
scattering probabilities W shows that for all reasonable types of inhomogeneities the decay of the surface-driven 
scattering probabilities Wjj> at large j,f is a function of (j + j') 2 R 2 / L 2 and (j — j') 2 R 2 / L 2 . Therefore, for small 
inhomogeneities, R/L <C 1, the the mode-coupling scattering probabilities Wjf with j' ^ j have the same order 
of magnitude as for the intraband scattering Wjj. What is more, the contributions of the higher modes to particle 
transport are quite noticeable and decrease rather slowly with increasing This is illustrated in Figure 5 for the 
power-law inhomogeneities (R/L — 0.1) where three curves correspond to diffusion in the single-mode, two-mode, and 
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three- mode regimes including the mode coupling. The surface correlator is exponential in momentum space, Eq. (|19|l . 
with /x = 1 (the Staras correlator). All three curves are normalized by the single value of <i(50) for the single- mode 
curve. As one can see from the Figure, the turning on of the mode coupling does not lead to any qualitative changes 
and results simply in the increase of the overall scattering cross-section. The contributions from all modes have the 
same order of magnitude. 

The situation changes dramatically at large R/L. Figure 6 presents the same calculation as for Figure 5 but for 
R/L = 100. As one can see, at small x = L/A all three curves coincide. This means that the contribution of 
higher modes and the mode-coupling effects are negligible even though the interband transitions are allowed. This is 
explained by the fact that in high-quality films with R/L ^ 1 the mode-coupling transitions switch on one by one 
when the values of the film thickness reaches the values (??) (with a small logarithmic correction which depends on 
the correlation function of inhomogeneities; Eq. (??) assumes that the quantum numbers j are small and qj ~ ft/A). 
Therefore, at thicknesses L < L\ (the first branching point in the curves) all interband transitions j <-» j' are 
suppressed. Since at R/L ^ 1 the absolute values of (jR/L) 2 grow very rapidly with j, the contributions of the 
higher modes are negligible as well, and the system is effectively in a single- mode regime. At L ~ L\, Eq. (??), the 
mode-coupling probability Ww becomes comparable to W±\ and the mode-coupling between the two lowest modes 
becomes robust. Only at this point contribution of the second mode become noticeable. Therefore, at L\ < L < L 2 
the system is effectively in a two-mode regime, at Li < L < Li - in a three-mode regime, and so on, as it is seen 
clearly in Figure 6. 

The results for the Gaussian correlator are similar to those for the power-law one. The exponential correlators, on 
the other hand, do not exhibit behavior similar to Figure 6. The power spectrum for such correlators is decaying very 
slowly at large q thus ensuring robust mode-coupling transitions at any R/L. As a result, the curves d(L) and a (L) 
resemble those on Figure 5 at any value of R/L. 

The important conclusion here is that in high-quality films R/L 3> 1 with the surface-driven scattering the con- 
tribution of the higher modes becomes important only as a result of the emergence of mode coupling at the values 
of the film thickness (??). However, after the mode coupling turns on at certain values of the film thickness, the 
contribution of the higher mode becomes much bigger than it is usually believed. 



E. Particles with non-quadratic energy spectrum 



It is very interesting to check to what extent our analysis of the mode coupling effects is sensitive to the form 
of particle spectra. We start from deriving an analog of Eq. (??) for non-parabolic particles, i.e, from finding the 
critical values of thickness at which the mode coupling become noticeable for individual modes in high-quality films, 
R/L > 1. 

Let us consider particles with a bulk spectrum e(p), or, after quantization, ej (q) = e (jrhj /L.q). Scattering by 
inhomogeneities of the lateral size R change the lateral momentum by Sq ~ h/ R. This small change in momentum 
(R ^> L) is sufficient for the interband transition j j + 1 only if the energy conservation, ej (q) = e J+ i (q — fr/R), 
can be satisfied: 

The solution of this equation gives the values of the critical thickness Lj at which the mode coupling channels j <->■ j +1 
become open. [More accurate equation than (|36[) contains coefficients of the order of 1 which depend on detail of the 
correlation function of inhomogeneities]. In the case of parabolic inhomogeneities, as for all particles for which the 
energy spectrum e (p) depends only on the absolute value of momentum p, Eq. (|36[) reduces to (??). 
As an example, we consider " ultrarelativistic" particles, e = cp, i.e., 

e 3 (q) = cy/frhj/Lf + q*. (37) 

The applications include photons between two (rough) mirrors or quantized phonons in helium films. For non-parabolic 
spectra e (p) , the equation for the scattering probabilities Wjj< should be modified as [4| 



(q, q') S (e m - e iv ) = ^ ^ - ^ + C 22 + 2( 12 (-l) i+/ ] S (e iq - e jv ) , (38) 



while qj/m in the equations for the single-particle diffusion (|12|l should be replaced by the lateral velocity vj — dej/dq. 
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Then straightforward algebra leads to the following expression to the single-particle diffusion coefficient: 

ci? 3 

D = —d{uL/Kc), (39) 

where Tjj> 1 is the dimensionless inverse relaxation time matrix @ with the changes mentioned above. 

Figure 7 presents normalized diffusion coefficient D (x) /D(50), Eq. (|39|l . for random inhomogeneities with the 
Gaussian correlation function; x — luL/ttc is the dimensionless energy (frequency). The curves are marked by the val- 
ues of R/L. The curves with small R/L exhibit the "standard" saw-tooth QSE. The curves for larger inhomogeneities 
exhibit the large-scale oscillations with the peaks which correspond to opening of the mode-coupling channels and 
which are described approximately by Eq. (|36|) . The curves for the exponential in momentum space surface correlators 
(|19|l are, essentially, the same. If the surface correlator has a power-law shape in momentum space, which means that 
there exist inhomogeneities of all sizes, the transport coefficients for particles with non-quadratic spectra assume the 
same " standard" saw-tooth shape as for the particles with a parabolic spectrum. 

The general conclusion is that the mode coupling effects, which manifest themselves in the new type of QSE in 
transport in high-quality films, are very robust and are not sensitive to the type of the energy spectrum. This means 
the this type of QSE should exist for (quasi-)particles of different nature and for various solid state systems. 



F. Interwall interference effects 



What also makes the scattering by surface inhomogeneities different from scattering by bulk fluctuations or impu- 
rities is the possible interference of particle reflected from the opposite wall. This interference is especially interesting 
in the case when the inhomogeneities from the opposite walls are correlated. In this Section we discuss the effect of 
this interwall correlation of inhomogeneities on mode coupling. The existence of this non-trivial effect is a unique 
feature of surface scattering that does not have any analog in scattering by bulk inhomogeneities. Surprisingly, the 
possibility of cross-wall correlation of surface inhomogeneities from opposite walls gives an interesting insight into 
mode coupling. The study of the effect of interwall correlation of inhomogeneities has been initiated in Refs. [a. l2f]l l2l| 
(for additional results in application to excitons see Ref. |3fJ). 

The effect of interwall correlations has two unique features stemming from the sign of the interference of scattering 
from opposite walls. Because of the (— iy +J factor in the interwall contribution to the scattering probability, Eci. (|23[l . 
the contribution of the interwall correlation function £ 12 has different signs for in-band (j — j') and mode-coupling 
(j = j' ± 1) scattering processes. Depending on the magnitude of £ 12 and its sign, its contribution can enhance or 
suppress the mode coupling effects. 

To decrease the number of parameters, we assume that, as in Refs. HIH EI, the correlation functions of inhomo- 
geneities on both walls Cn an d C22 are gi ven by the same function, ( u (s) — C22 ( s ) = C( s )- The structure of the 
interwall correlator of inhomogeneities ( 12 (s) is assumed to be the same as for the intrawall correlations with the 
same correlation radius R. However, the amplitude a of the interwall correlations is different from the intrawall ones, 



C11 = C22 - C (s) , C12 (a) = «C 00 , |o| < 1. (40) 

Note that in contrast to the on- wall correlation functions C11 •> C22 : t ne sign of the interwall correlation function 
£ 12 is not fixed; even £ 12 (s = 0) can be negative. By itself, the sign of the interwall correlations £ 12 = • £ 2 ) 
is ambiguous and depends on how do we introduce the signs of the deviations of the wall positions ^ 2 from the 
averages ±L/2; throughout this paper, we use the definitions l|14|) . With this definition, the sign of a can be positive 
or negative depending on whether the inhomogeneities from the opposite walls "attract" or "repel" each other. If 
the inhomogeneities from the opposite walls simply reproduce each other ("parallel" walls; the film thickness is 
constant along the film), then, with our definition of the wall inhomogeneities Ea. i|14|) . a = — 1. This type of interwall 
correlation is likely to occur when an ultrathin film grows on an inhomogcneous substrate. In the opposite case of 
walls with opposite modulations (" antiparallcl" walls), a = 1. This is the case of pure thickness fluctuations, which is 
likely to occur, for example, after the film (wire) has been inhomogeneously stretched. In the case of "parallel" walls, 
a — > — 1, the destructive interference of scattering by opposite walls can, in the absence of mode coupling (see below), 
completely negate all transport manifestations of the wall corrugation. These two limiting situations is presented in 
Figure 8. In general, — 1 < a < 1. 
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To extract the effect of interwall correlations, we will present the results for the transport coefficient at different 
values of the interwall amplitude a and calculate the relative change of conductivity a and diffusion coefficient D 
caused by the introduction of such correlations, 

tip) - f (0) //(<») _ (f(0) 

<*(<») = L J — - (41) 

where and (d<-^ and d^) are the values of a and £>, Eqs. (|29|l , H30|l . calculated with and without interwall 
correlations. An additional benefit of this definition is that the functions <jr a ' are automatically normalized. 

In the presence of interwall correlations 14U|) . the transition probabilities Wjy (q, q') l|23l) become proportional to 

l + «(-l) J+/ ]c(|q,-q>|). (42) 

The most interesting effects of the interwall correlations are related to the oscillating structure of the term with a in 
Eq. 62J. 

If the mode coupling is suppressed, then the only important terms in Eq. (|42|) are the diagonal ones with j = j' 
and the function tfi^ is a constant, 

<p^ [no mode coupling] = 1 = . (43) 

1 + a 1 + a 

In this case the presence of interwall correlations leads to a simple increase or decrease, depending on the sign of a, 
of the transport coefficients by the factor — a/(l + a). This is always the case, for example, when only one mode is 
important [8(ij . Therefore, all deviations of (f>^ from the constant 1|42|1 are due solely to the scattering-driven mode 
coupling. This gives a non-trivial insight into the mode coupling and its consequences. 

For example, as it is clear from Figure 6, the mode coupling for the power-law inhomogeneities with R/L = 100 
appears only at x = L/A > 15. Therefore, the function (j>^ (x < 15) should be flat, Eq. (|43|l . and exhibit anomalies 
in the points in which the mode coupling effects are switched on. This is illustrated in Figure 9 which presents 
the function <f>^ a ' (x) for the same power-law inhomogeneities as Figure 4 (R/L = 100) at five different values of a, 
a = —0.5; —0.1; 0.1; 0.5; 09. The flat parts on left hand side of all curves (x < 15) correspond to the absence of mode 
coupling at these values of R/L and are given by Eq. I)43[l. The peaks in the curves show the values of x = L/A 
for the consecutive openings of the mode scattering channels. The difference in the amplitudes of the peaks is easily 
explained by the dependence of the scattering probabilities on the interwall correlation amplitude a, Eq. (|42|l . 

The contribution of the term with a in l|42(l has a different sign for different Wjf depending on whether j + j' 
is even or odd. Since the mode coupling channels j < — > j + 1 in high-quality films with Gaussian and exponential 
power functions turn on one by one with increasing L, one would expect that the function <fi^ (L) in such films 
should become a step-wise function of L. This is not correct. In films with large R/L and frozen out mode coupling 
effects the contributions of individual modes decrease rapidly as 1/j 4 [2(j. However, when the transitions j < — > j + 1 
are switched on, the overall contribution of the mode j + 1 increases disproportionately (Figure 6). As a result, the 
function tp^ becomes an oscillating rather than step- wise function as it is seen clearly in Figure 9. 

The positions of the QSE peaks in systems with interwall correlations differ from Eq. I|35|) and depend on the 
value of a. The shifts of peaks in Figure 9, which depend on the value of a, are better illustrated in Figure 10 in 
which we presented the normalized diffusion coefficient itself (and not the function 4>^) for the same type of surface 
inhomogeneities and the same value of R/L = 100 for two different interwall amplitudes, a = —0.9; 0.9 as a function 
of x = L/A. The explanation of these shifts is the following. Let us assume that the first peak is observed at x = x\. 
In this point the value of W\2 (x) reaches Wu (x), W\i (x\) — W12 (xi) . According to Eq. [fl^ , in the presence of 
interwall correlations, these scattering probabilities change with respect to their values in the absence of 

interwall correlations as 

W[f = (l + a) , W[f = (1 - a) W[f . (44) 

Since near the peak position (x) grows very rapidly while Wyy does not change much, ( x ) reaches the 

value wffl earlier than (x) reaches at negative a and later at positive a. This exactly what is happening 

in Figure 10. At small a, the change in position x± with respect to its value in the absence of interwall correlations is 

Axi = 2aW^ (asi) / \dW$/dx ~ dW^/dxl . (45) 

The oscillating nature of the interwall contribution, Eq. I|42|) should be more pronounced for the systems with 
smaller inhomogeneities, in which the mode-coupling transitions are as probable as the intraband scattering. In this 
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case the flat areas (|43f) should be absent. Instead, the curves <jy- a ' (x) should exhibit QSE with the oscillations in points 
xs in which the number of occupied minibands 5 changes by 1, 5 — > 5+1. [With our choice of dimensionless variables, 
the period of these oscillations is equal to 1]. This is illustrated in Figures 11 and 12 which contain the data similar to 
those in Figure 9 but for smaller inhomogeneities, R/L — 0.1 and R/L = 1 respectively. At R/L <C 1, when all mode 
coupling transitions and intraband scattering are equally probable, the addition of an extra band 5 adds all interwall 
terms with the sign-changing coefficients a (— lY +S , Eq. i|42|) . However, since the main mode, j = 1, contributes the 
most to the transport flow, the overall sign of the interwall contribution is the sign of a (— 1) 1+S and should change 
in the points in which 5 (x) changes. This is exactly what can be seen in Figure 11. The amplitude of the oscillations 
grows with an increase in \a\ and goes down with increasing x. Figure 11 (R/L = 0.1) demonstrates these oscillations 
for large interwall correlations, a = —0.9; 0.9. Since the signs of these two interwall amplitudes are opposite, the 
contributions from these two types of cross-correlations have opposite signs, Eq. (|42|l . The analogous curves for all 
interwall amplitudes \a\ < 0.9 are squeezed between the curves for a = —0.9; 0.9. At larger inhomogeneities, R/L = 1 
(Figure 12), one can still see the well-pronounced QSE oscillations, but the average is already noticeably shifted from 
zero as it should be at larger R/L (cf. Figure 9). 

Note, that the height of the first peak is always given by Eq. ()43[) and can be quite large when a — > — 1 ("parallel" 
walls). At a — » — 1 the interwall correlation compensates almost completely for dephasing caused by scattering 
from individual wall inhomogeneities. In this case, if the wall scattering is the only relaxation mechanism (ballistic 
transport), the lateral mean free path becomes infinitely large even if both walls are rough! In Figures 11 and 12 the 
height of the first peak for a = —0.9 is 9 and the peak does not fit into the Figures. For positive values of the interwall 
amplitude a (" antiparallel" walls), the second peak has the largest amplitude, while the first (negative) one, which is 
given by Eq. (|43l) . has a smaller amplitude. 

In conclusion, the possible correlation of random inhomogeneities from the opposite walls provides a non-trivial 
insight into the mode coupling. On the other hand, measurements of the dependence of the transport coefficients on 
the film thickness or particle energy can provide unique information on the interwall correlations since, depending on 
the situation, the effect of interwall correlations can be constructive, destructive, or oscillating. The shift of oscillations 
gives the information on both the strength and sign of the interwall correlations. 

G. Scattering by multiscale inhomogeneities 

Above we studied the systems with random inhomogeneities of a single, well-defined spatial scale (correlation radius) 
R. In the case of single-scale inhomogeneities, such as inhomogeneities with a Gaussian or exponential power spectrum, 
the mode coupling channels at large R open one by one at definite values of the film thickness (|35[1 leading to large scale 
QSE oscillations of the transport coefficients. In the opposite case of the power spectrum with inhomogeneities of all 
sizes, such as slowly decaying power-law power spectrum with a low index, the mode coupling is always robust leading 
to the disappearance of the large-scale oscillations and the restoration of the standard saw-tooth QSE. It is interesting 
to investigate the behavior of QSE in an intermediate situation in system with inhomogeneities of few distinct scales. 
Figure 13 presents the data for the diffusion coefficient for a film with Gaussian surface inhomogeneities of three types: 
the inhomogeneities with R/L = 25, R/L = 10, and the inhomogeneities with a combination of both sizes (the sum 
of the corresponding power spectra). All three curves are normalized by the same value d {R/L = 25; x = 50) and are 
labeled by the value of R/L. As one can easily see, the combining the inhomogeneities with two correlation sizes does 
not lead to a mechanical mixture of the individual oscillations but results in smoothing, shifting, and rescaling of the 
oscillations. Adding several more scales leads simply to a disappearance of the QSE oscillations. At present it is not 
clear yet the combination of how many scales is necessary for the restoration of the saw-tooth behavior. 

VI. CONCLUSIONS 

In summary, we analyzed the mode coupling and its effect on QSE in transport in high quality quasi-2Z? quantum 
systems (R/L ^> 1) with various types of surface, thickness, and bulk fluctuations. Here are the conclusions: 

• Though the transport equations and mode coupling effects for systems with bulk and interface fluctuations look 
almost identical, QSE in such systems is not the same. The appearance of large scale oscillations of the transport 
coefficients requires not only the opening of mode coupling channels at distinct values of film thickness (I35|l . 
(|36|l as it happens in both surface- and bulk-driven systems, but also the predominant role of modes with low 
quantum numbers (grazing particles) . The latter requirement is routine for surface scattering but is not fulfilled 
for bulk fluctuations of a general form. Only if the fluctuations in the bulk do not depend on the coordinate 
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across the film, Eq. (|25fl . the transport coefficients manifest the same large-scale oscillations as in the case of 
surface scattering. This explains a huge difference in QSE in high-quality films with bulk and surface scattering. 

• One of the most striking conclusions concerns the contributions from different modes in high quality samples in 
the case of surface scattering by large inhomogeneities, R/L ^> 1. Without mode coupling, the contributions from 
a mode with quantum number j would be proportional po| to 1/j 4 and higher modes would have been almost 
irrelevant. These higher modes contribute to transport only because of the mode coupling. This conclusion is 
of little interest for small-size defects R/L <C 1 since in such systems the scattering-driven mode coupling is 
always robust and the contribution from the higher modes is important. 

In high-quality films R/L 3> 1, the mode-coupling transitions switch on one by one at the values of the film 
thickness (I35|) . (|36|l . Thus, the higher modes become important also one by one, only after the corresponding 
mode coupling channel is turned on, Figure 6. After the mode coupling turns on, the contribution of the 
higher modes is much higher than one usually assumes and the description that singles out the grazing particles 
becomes wrong. 

• The consecutive opening of the mode coupling channels in high-quality films at distinct values of the film 
thickness, Eqs. (|35l) . (|36[) . which leads to giant QSE oscillations, is a very robust effect that is not very sensitive 
to the nature of (quasi-)particles and the form of their spectrum. As a result, the effect can be observed in a 
wide variety of quantized systems such as metal or semiconductor films, quantum wires, ultranarrow channels, 
optical fibers, etc. 

• An interesting manifestation of the coupling effects in high quality films is related to possible correlation between 
random inhomogeneities from the opposite walls. The interference of scattering from the opposite walls changes 
its sign from constructive to destructive depending on the parity of the sum of mode quantum numbers and, 
therefore, provides contributions of the opposite signs for intramode and mode coupling channels. In some cases 
("parallel" walls) the opening of the mode coupling channel can be responsible tor the cutoff for the mean free 
path for grazing particles which would be nearly divergent otherwise. The presence of interwall correlations 
can help to distinguish films with surface and thickness fluctuations. The shift of conductivity or diffusion 
oscillations provides the information on both the sign and strength of interwall correlations. 

• The presence of multiscale inhomogeneities with several distinct correlation radii R leads, instead of a mechanical 
mixture of individual QSE patterns, to shifting and smoothing of the QSE oscillations of the transport coefficients 
that are inherent to high-quality films with a single-scale roughness. 

• The results can lead to new, non-destructive ways of studying the quality of the high quality surfaces, including 
the buried surfaces and interfaces, by measuring the lateral conductivity or diffusion. This is especially valuable 
for high-quality surfaces with large-scale inhomogeneities for which the usual scanning techniques can become 
problematic because of very large scanning areas. 

The results of the paper can be applied to particles in a wide range of quantum quasi-2D systems. The results 
can also be cautiously extended even to quasi-l-D systems up to the point when the strong localization effects render 
transport calculations meaningless (see review [3ll | and references therein). 
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VIII. FIGURE CAPTIONS 

Figure 1. (Color online) Normalized conductivity / (x) // (50) , x = L/Ap, Eq. I|29[l .with artificially frozen mode 
coupling. The shapes of all five curves are identical irrespective of the type of inhomogeneities (Gaussian, exponential, 
or power-law) and the value of R/L. With this resolution all the curves are smooth. 

Figure 2. (Color online) Normalized single-particle diffusion coefficient d(x) /d(50), x — L/A with artificially 
frozen mode coupling. The shapes of all five curves are identical irrespective of the type of inhomogeneities (Gaussian, 
exponential, or power-law) and the value of R/L. Small saw-tooth anomalies correspond to changes in the number of 
accessible mini-bands S(L), Eq. 

Figure 3. (Color online) Normalized single-particle diffusion coefficient with surface and bulk scattering; in both 
cases the inhomogeneities are Gaussian with R/L = 50. 

Figure 4. (Color online) Normalized single-particle diffusion coefficient with surface and bulk scattering; in both 
cases the inhomogeneities are Gaussian with R/L = 0.1; 1; 10; 50. For surface scattering, all four curves, marked by 
the values of R/L, are different. For bulk scattering, there are two groups of coinciding curves with small and large 
R/L. 

Figure 5. (Color online) Normalized single-particle diffusion coefficient d(x) /d(36) , x = L/A for power-law inho- 
mogeneities, Eq. I|19|) at fi = 1 for small size inhomogeneities, R/L = 0.1. Curve 1 takes into account only the main 
mode. Curve 2 accounts for the first two modes, including coupling, curve 3 - the first three modes. It is clear that 
all three modes are equally important. All three curves use the same normalization parameter d(36) taken from the 
single (main) mode contribution (curve 1). 

Figure 6. (Color online) The same as Figure 5 but for large-scale inhomogeneities, R/L = 100. The splits occur in 
the points when the mode coupling becomes noticeable. It is clear that the contributions of the higher modes become 
noticeable only when their coupling to the main mode becomes large. 

Figure 7. (Color online) Normalized diffusion coefficient D(x)/D (50), Eq. (|39|l . for ultrarelativistic particles, 
e = cp, and random inhomogeneities with the Gaussian correlation function. The curves are marked by the values of 
R/L; x — ujL/ttc is the dimensionless energy (frequency). 

Figure 8. (Color online) Two different types of walls with correlated random inhomogeneities. For "parallel" walls 
the interwall correlation amplitude a = — 1, for " antiparallel" walls - a = 1. In general, the interwall correlation 
amplitude — 1 < a < 1. 

Figure 9. (Color online) Relative contribution of the interwall correlations to the single-particle diffusion, Eq. ()41JI . 
Large-size (R/L = 100) power-law inhomogeneities, Eq. (|19l) with /i — 1. The curves are labelled by the values of the 
interwall amplitude a. 

Figure 10. (Color online) Normalized single-particle diffusion coefficient d(x) /d(36) for the same inhomogeneities 
as in Figure 8 for two values of the interwall amplitude a, a = —0.9; 0.9. 

Figure 11. (Color online) The same as in Figure 8 but for small inhomogeneities, R/L = 0.1. The curves are labelled 
by the values of the interwall amplitude a = —0.9; 0.9. 

Figure 12. (Color online) The same as in Figure 8 but with R/L = 1. The curves are labelled by the values of the 
interwall amplitude a. 

Figure 13. (Color online) Diffusion coefficient for Gaussian inhomogeneities with R/L = 10, R/L = 25 and with 
the sum of inhomogeneities of both sizes. The curves are labled by the value of R/L. All three curves are normalized 
by the value of d (R/L = 25; x = 50) . 
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